METHOD AND APPARATUS FOR NOISE FLOOR ESTIMATION 



Field of the Invention 

The invention relates to the field of telecornmuiiications and in particular the signal 
processing in receivers of cornmunic&tion signals. 

Background 

In the art of telecommunications, it is desirable to estimate the average noise power 
contained in a received signal. Typically, the received signal is a frequency division 
multiplexed composite signal containing several distinct signals. Jt is important to be able to 
estimate the noise power without a priori knowledge of the number, nature, carrier-frequency, 
bandwidth or amplitudes of the distinct signals. 

In the art* noise is considered an additive phenomenon. It distorts a signal by adding 
random time- varying values to the amplitude qf the non-noise component of the signal. It is 
thus of great importance to develop a radio receiver that is robust to the noise phenomenon. 

The general problem of estimating the variance of the additive noise in a communication 
signal composed of many narrowband frequency division multiplex signals is important in 
many radio communication systems- It allows setting the detectors threshold to its maximum 
sensitivity while limiting the probability of a false alarm to a value less than or equal to a 
desired leveL The detector decides that a signal is present when the energy of a signal, over a 
period of time, exceeds a threshold value. The detected signals can then he passed to a 
receiver or to measurement equipment, for demodulation or to acquirer more information 
about the signal. 

Currently there are several techniques that are used in this pursuit, but all come with 
sacrifices. One such technique is to reserve a portion of the bandwidth of the composite 
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signal to cany a "noise channel" Thus by measuring the power in this reserved band, an 
estimate of the noise power in the other channels can be obtained. It is common though for 
these "noise channels" to be difficult to identify as their placement in the spectrum of the 
composite signal is not governed by any global convention which leads to different channel 
assignments based on geographical location of the measurement equipment Other errors are 
introduced when the "noise channel" is affected by the signals in adjacent channels. This 
results in a biased estimation of the noise floor. If the channel assignment is unknown,, it may 
not be possible to identify the 4 'noise channel" and determine the noise floor power level. In 
addition to these problems the noise floor tends to fluctuate with environmental and operation 
co nditi ons. 

In the art, it is common practice to set a noise threshold value, determined in large part by 
the estimated noise floor, whereby a signal is considered present if and only if it is above the 
noise threshold. 

The traditional approach is to set a threshoLd value depending on the user needs. When 
high sensitivity is required, the threshold is set to a value close to the user observed noise_. 
floor leve l. When only strong signaLs are of interest, the threshold is set to a high value 
relative to the observed noise floor leveL This approach is highly subjective and may result in 
being too aggressive or too conservative for the actually measured signal environment . The 
threshold is typically set manually in a communications system. This is not feasible when 
fast response to dynamically changing conditions is required. Very few automated 
techniques have been studied to address the general problem of rapidly and dynamically 
estimating the noi se floor level or the si^lJo^noise.ratio. One typical method is to isolate a 
channel and to assign this channel as a noise only channel for estimation purposes. This 
method is difficult to apply in many new digital communication systems where the resources 
are often utilized close to their capacity. It is also a problem for dynamic signal environments 
based on multiple access. Another typical method is tq assume that the composite wideband 
signal non-noise subspace dimension is known. This is again difficult to realize in practical 
situations with unknown signals where robustness is a key feature. On another front, a 
technique based on morphological binary image processing operators (similar to rank-order 
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filtering) as described in "Automatic Noise Floor Spectrum Estimation in the Presence of 
Signals", MJ. Ready, M.L. Downey, and L.J. Corbalis, Proceedings of the 31" Asilomar 
Conference on Signals, Systems and Computers, vol. 1, pp 877-881, November 1997, 
processes a binary image of the wideband received power spectrum of a composite signal has 
been proposed. Thus, the approach does not process the time samples or frequency domain 
transformed data directly, but the instead it processes an image of the spectrum. This involves 
two dimensional signal processing that is usually complex. Furthermore, practical 
performance results have yet to be quantified with such an approach. In general, techniques 
based on textbook results tend to assume several conditions that are not truly present in many 
bands of the radio spectrum. 

The computational complexity of a method is a limiting factor of the method's 
applicability. Typically, a spectrum monitoring receiver cannot observe the full spectrum it 
is allotted, thus it divides the spectrum into smaller blocks that can be scanned, each of these 
blocks possibly contains a plurality of signals. Each block is received sequentially, an 
estimate of the noise floor is arrived at, and then the receiver shifts its frequency to receive 
another block. This cycling through the allotted bandwidth must be done rapidly to ensure 
that all messages being transmitted in the different parts of the allotted spectrum can be 
received. If the method to derive a noise floor estimate is too complex then it either cannot 
keep pace with the shifting of the block or it causes the shifting of the blocks to slow down, 
in the event of either result the method ceases to be useful. 

Thus, it would be desirable to design a system that resolves the problems that the above 
mentioned methods have with narrowband signals of unknown nature, as described in "Noise 
Floor Estimation for Wide-Band FFT Filter Banks". It would also be desirable for a new 
system to have fully digital characteristics so as to permit real-time signal performance, 
without requiring the dedication of a reserved frequency band for use as a noise channel. 




Summary of the Invention 

It is an o"bjcct of the invention to remedy numerous problems found in the prior art 
methods. 

It is another object of the invention to provide a noise floor estimation method and 
apparatus capable of operation on a wideband signal without full knowledge of the nature of 
the narrowband signals contained therein. 

It is another object of the invention to provide digital characteristics to allow real-time 
signal processing with low levels of computational complexity. 

The preferred embodiment of the method of the invention, with very little a priori 
knowledge of the signals, allows both robust performance in non-ideal environments, and an 
automated maximum sensitivity threshold setting method. These modifications make the 
process of the invention an attractive solution to problems emerging in the field of digital 
communications. As an added benefit, information derived by the method of the invention 
may be used in spectrum quality indication, interference analysis, and signal demodulation. 

In the preferred embodiment of the method in accordance with the present invention the 
noise floor is estimated by carrying our the steps of: 

a. representing the continuous signal as a series of discrete frequency and 
amplitude values; 



b. creating a histogram based on the discrete frequency and amplitude values; 
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c. deriving a no ise floor estimate from characteristics of the histogram. 
The discrete frequency and amplitude values are preferably obtained hy the steps of: 

d. sampling a received wideband signal by a plurality of analogue-to-digital 
converters to generate a series of output signals; 

e. windowing the output signals of the analogue-to-digital converters; 

f. applying a mathematical transform, preferably a Fast Fourier Transform, to the 
results of the windowing to obtain a series of discrete frequency values; 
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converting an amplitude of each discrete frequency value to log-domain 
representation; and 



h. rounding the log-domain representation of the amplitude for each discrete 
frequency value to the nearest integer value to generate a discrete amplitude 
value. 

The windowing step preferably includes the steps of: 

i. selecting a discrete weighting function; 

j . multiplying the value of each output signal of the series by a corresponding 
element of the discrete weighting function. 

The amplitude ofeach value of the discrete frequency series is preferably converted to log 
domain representation by multiplying 20 by the base 10 logarithm of the magnitude of the 
element; whereby the log domain representation of the amplitudes results in the amplitudes 
being expressed either as decibel (dB) values, or as decibel milliwatt (dBm) values. 
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in a first preferred embodiment of the steps of crying and analyzing the histogram to 
eshmate the noise floor, herein referred to as the discrete level crossing me method, a 
histogram chart is constructed where the minimum value is equal to the minimum mteger 
value in dB of the received signal strength, and the maximum value is equal to one less than 
the highest integer value in dB of the received signal strength. Every time a segment of the 
signal is found to have a positive slope, the histogram bins crossed by the segment are 
i—ted by on. Thus, for example, if a given segment starts at -^^"^ 
dB in menextsample, the values ofth. bins representing -11®, -1,9, -1**, -107 and -1,6 dB 
are incremented by one. After the complete Mstogram has been assembled, the next step ,s 

) p^blytodete^^^ 

histogram will be bimodal in the presence of noise and aplurality of signals, wvm the first 
true local maximum indicating the noise floor. To find the noise floor, preferably the first 
local maximum string from the minimum quantised dB value is selected. Because me 
. histogram is created from sampled data rtbetedby a random source it is possible to have . 
5 numerous !ocal maxima that do not reflect the true noise floor, as a result it * preferred to 
check if the local maximum is in fact the n^floor. In the most preferred embodiment, the 
method for the maximum value found includes the step of checking if a histogram value 
within Y dB is higher than the current maximum, where Y is a preferred number denved 
empirically from the characteristics of the system, most preferably between 2 and 3 dB. If 
20 the no value in the next Y db is greater than the value at the maximum then the current 
maximum abscissais indicative of the noise floor level estimate per bins, however xf Hie 
value is greater Y db away from me maximum than it is at the maximum then it is clear that 
the noise floor has not been found. In the event that the maximum is not the noise floor, the 
ne xt local maximum is preferably used and theYdB test ap P Ued again until the test faUs to 
25 produce a value higher th*x the value at the maximum. The obtained abscissa value is finally 
preferably corrected by a fixed offset ofXdB, where X ,s preferably dependent on the 
selected window, to give the noise floor level estimate per bin for the processed block of data. 

In a second preferred embodiment of the steps of creating and analyzing the histogram to 
estimate the noise fl.«r, herein referred to as the penalized quasi log^elihood method, a 
30 histogram of quantized frequency domain amplitudes is constructed with each but 1 dB m 
size The histogram is then transformed mto a sorted vector of linear amplitude elements m 



decreasing order. The sorted imea, vector is P« *° m *» ™ to 

rZl S a^po lK co„ 5 ecudvess r o ple sof t hesor tt dvec.o,I n .h I s I n el hod M , S 

I " hTarglmanr, ^ log-likelihood function L<k) i. *-^*-*r*« 

5 pref«rab,y added ,„ . penaity prJy— function pW **• - =F~» ° f "» ^ 
-L« + C „(*). ^ere C is a constant, preferably empirically deurottfed and most 
prerer^vnear - 2 .o. TUe m — of Ute-ZW^pW^cdoni, preferably =e,ee K d 

wiererheWextoftheminirn^amplimdeisdenotedby U a line is drawn on the 
Histogram..^ ,„,itca»be seen that the received signal be divide. tn to2 
10 groups, the first betas the menage subspace, wherethe informal contained ^.heoon- 
„1 sign., exists, and *. noise subsp.ce, where the random n,,e 
noise «£ level estimate per frequency bta is preferably obtalnedoy divtdtng by K the 

M . _ismdtetval„e S of^Ms<«=dgroup.Thu S the™isefl«>restir„ M e P er 
average M - q SF - Lsm^ucat 
^ K^^^ - the average value of the power along the pensions of ncse subspace. 
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Th.appueatiouwheremcproposedmeu^^ 

^ns are necessary and - adaptive thread v*ue is quired. The, can also be u~d to 
1— .engthymanua, noise .eve! fl< "" 
level estimates can be used as an indication of the spectrum quality. 

As stated earlier, it is anobject of the invennon to derive an estimate in real-toe of the 

terms of carrier «y, baudwidm. and power, and using only . shon-nme s,gna>. feft. 
preferred apparatus in accordance with the invention, the wideband signal, —tg a 

7 „ i™fer*lyreceivedfromtheanalogueoutp»tofarad,o 

olirrahty of narrowband signals, is pieteraoiyrev 

The sampled amp.ih.des, though discreK, are real valued and thus the s.gn^ ts sUU 
considered OOBlllIM . ta ^^^^^^^^ILl^. 

collection of ccmp.e* ampHtudes at several discrete feq uencies. The amphtudes of 



signals at the discrete frequencies arc converted in a converter module to the log domam 
representation (decibels (dB) or decibel milliwatts (dBm)) «nd rounded to the nearest mteger 
-values. Thus the received signal is now represented by bo* discrete amplitude and 
frequency valves. 

5 Brief Description of the Drawings 

Exemplary embodiments of the invention will now be forth* described with 
references to the drawings in which same reference numerals designate similar parts 
throughout the figures thereof, and wherein: 

Figure 1 is a block diagram of the signal processing elements for the discrete level 
10 crossing rite (DLCR) embodiment of the method according to the invention; 

Figure 2 is an example of the discrete level crossing rate with a noise only input 

signal; 

Figure 3 is an example of the discrete level crossing rate with both noise andnon- 

noise signals at the receiver's input 

1 5 Figure 4 is an example of the behaviour of the DLCR method as a function of the 

occupancy; 

Figure 5 is an example of experimental results showing the time behaviour with low 
and high signal frequency band occupancy for the DLCR method in the 800 MHz band; 

Figure 6 is a block diagram of the signal processing elements for the penalized quasi 
20 log-likelihood (PLLM) embodiment of the method according to the invention; 
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Figure 7 is . Plot of the ouasi I***- ««*- «" — ^ """"^ 

polynomial; 

is an example of the PLLM(k) Auction for anosie only signal, and a nan- 



Figure 8 
noise with noise signal; 

Figure 9 is - example of .hebehaviour of 0»PLLM ro ethod a, a fhncSon of the 

frequency band occupancy; 

Figm = io is an of experimental —* * e Hme b *" h " " i,h 
ar»> high sign* occupancy for the PLLM method in the 800 MB* band; and 

Table 1 shows the offset for the DLCR method. 
,0 Detailed Description of the Preferred Embodiment 

!n the context <*« 
are observed by a wide-band ehamteUsauon receiver. The receiver rentes 

when Ore only assumption is an input signal with additive noise nonnaHy debuted. For 
^ed informanoo vector, which is provided by an anaiogue-^K^ converter. 



samp 



T ta ( mw ^^»i^^*^ w "^ 4v, * ! ' ,,l ° ,I<h, , 

r.soiutionreauired, the signal processing stages *at follow and *. physical —ts of 
moving the data to the subsequent stage. To ensure a fas, scanning operauon of the 

be town by the monitoring equipment. The a ^ Knowledge of the recover . «n> 
and is limited to the information that signals are not wider than the samphng bandwdth. 
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Fqt a te - Wl— of « ntme * 0d , to - *e 

„ , *W thesi^lpathofthecstiinationrnelhQdispresentedmFxgurel.The 
Crossing Rat, method, the s lg nal patn ^ ^ ^ 

^.Ugita, — 10, ^ sign.. - stores . — of - - • 

V, i, ,h. vector r 103 The vector r 103 is then windowed or polyphase 

ri „ erl on specific requirements relive to . partite signal environment or or, SW 

signals)upt ,„ 0 ,h«taKsi Sm alpoworatthcailteDnainput.Iti5therofore 
of to filter bank . represent.^ the true «-l P» 
m^e^toeosiireftattheoquipinentispropetlyeaUDratefltoavotaany 

SZL floor .evei. » the dtgita. - « - - * " Til « 

oi uu« rnt5rtMfl i pnmrtlexitv but other transforms are 

, < r„, ir ; rt Transform (TFT) to minimize computational complexity, oui 

reprLntedby the vector* .05 is — * "■^■^"^ 

inlgeriog rnagmtude domain (dB or dBm) with a rounding operanon on20«nesthe 
io^baselOof^roagrut^ofeachcon^^pointofutev^The 

20 ^eroentattonofUte^— S ,e P 106 * simplified by the round.* ***** 

T^eiac* S ingRate(LCR)<unctionisweU^ _ . 

WCY Lee MeOraw-Hill, 1982. TheLCRis a s*=ond«>rder staasucal operants that ,s 

m the luation of me funcfio. The LCH isdefined as the tota, nurnber o^d 
30 clings of the signal envelope, a. a given amplitude A, per second. The standard LCR 

Z£ on time domain signals and is a continuous leve, function with a —us stgn* at 
II put. Tie current signal has a discrete-fluency discCe-anrpUn. de rep— on. To 
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c^TedlnThe upward direction. Therefore, wt^EM****** « 
rfle LCR is « - a his.ogram of the number of occurrences **- P-*» 
(positive slope) of. ^ amplitude vaiue in dB. For an additive Gaussian « only 
signal, this procedure produces a nmction havmg a shape of the form . 

V2pl0^ 
n[Z = A\= T"\ 



where,, is the ^er level in dB at which the number offing is recced, P*» 
correction parameter of the frequency leveis, and ^ 1. me va^ce of the additive complex 
Gaussian poise process in an FFT bin. The function above has a single ** can ^ 

«ny be verified to be located a. 10k*. (<r' ) . The parameter to estimate is the value of a 
and i, car.be estimated with the position of A. maximum of the function M\ when only 
noise is present 

The situation where one or several non-noise signs* are present with the additive 
„„ is e needs additional signal processing. For this case, it is possible that the maximum of the 
a^on is no, me noise floor level since me signal might cross . given ampUmde more often 

large bias errors, the fust local maximum, sorting from the lowest dB values, is selected. To 
chec* if this local maximum is no, due to . random Ducutatioain the constn.c,on of a ^ 
discrete LCR histogram, i, is then verified ft. there is no greater histogram value on the rrght 
side of the local maximum closer ,o Y dB relative to the curra* loed maximum- Thvs 
verification also stabilizes the method by preventing large ovemstima,e» of the nose floor 
level when the channel occupancy is high. The search for the appropriate maxrmum 
panned in search module 110, on me outpu, of the Discrete I*vel Crossing Rate module 
108. which is referred to as .09, « have different forms. The value - *" 

dynamic range of the non-noise signais. For mis apphcation, a value of 3 dB . used. Other 
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mC cha«s,ab U izetometod^^^^ 

possible The abscise v*,ue of to m— found in search module .10 ,s to notse floor 
eve, estimate 111 per FFT bin w«n an offset. The o t fse, .0 be added <otoabsc.se 
.iepends on 0- window o, to polyphase <«. used. Table , gives to offset for four poplar 
5 windows. Thus to no.se floor level estimate & 1 per FFT bin m dB .s 

3> [dB] = position of appropriate maximum [dBl + Offset [dB]. 
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The noise floor level estimate, <r car, be average for the am. frequency baud to 
reduce the variance if the noise floor level ts known to vary very slowly. 

To evaluate the performance of this method, simulation software was employed and a 
10 real-time implementation was deveioped to gamer static about the noise floor ieve! 
esHmates. The performance .» terms of the bias and the variance of the nc*e floor ievel 

sampled FM signal in the V39to 144 MHz band was used as me trauarmtted s.gna, m0.e^ 

baDi used a Blactonan window. The transmitted sisals all have the same av^ge power 

the receiver input. 

Figun: 2 demonstrate. to oonatntcuon of the DLCR histegran, for a particular vector 
x whenonry anoise signal is present in the simulation data at the input. InRgure 2«, a 

values^ resulting DLCR is presented in Figure 2h where Are nm^n Ml of the 
histogram is a. -* dB. Adding to 3 dB oftse, gives anoise floor level per FFT Mn of ^ dB^ 
to true value. Figure 3 demonstrates to —ion of to DIXRhistogram for a purhcular 
vector .when a_te having a channel chancy of 1..5V. 
i5 channel is preset with to noise signal in to sitnnlauon date a, to input of the A^The 
si^s have a signa. power 2 0 dB above to no.se floor revel. In Figure 3a, aBiaefcman 
^ i ^1024«tebrt^^»^^' U ^ Ml,B - 1,,t 
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* ■ pJcx \* where the maximum of the histogram is now at -8 dB 
resulting DLCR is presented in Fig. 3b where the maxi 

301 With the 3 dB offset, the noise floor level per FFT bin 

301 - . . TnEenCTa vth e algorithm will introduce a small D1 as as a function 

introduced in mis example- In general, the atgon 

^ a ,ienal spectrum shape, the signal amplitude and more 

nfihe narrowband non-noise signal spccuui v 

ot the narrowD foUowin(i the Dumber of channels M has been 

importantly the channel occupancy. In the following, rav 

ZT. Id *- — of bins p. channel K to 8. *~ — - . 

; siW8 J ~ p— « * ~. — ; c t * IoiK 

„ fl n ->ft m d 30 dB in 30 kHz with the true value of the noise 
values ofsifflial-to-noise ratios of 10, 20, ana ^uo»m j 

^ at 0 dB- Each column bar is the average of 1000 independent «* 
1 error bar for the standard deviation. Fisnre 4 shows 4a. *e esunrat* t— — i 

ceupancy 40! . Fisure 5 presents 1^1— — -*» — 

MHz bandwidth sampled at 5.12 MHz with 5 kHz bin re— The -*«™ 

,™.r,t,d - one where the centre frequency was at 848 MHz wrtn a 
thus 1024. Two curves are presented, one wnerem 

,„w occupancy of approximately 3 %. and one where the centre neoue^cyu 880.05 MB. 
IJe tjoccupaney is about 54 % with an averse s iS »al power of 30 dB above .he no.se 
wheretneoccy ,„„,„ each other it is assumed that the noise floor 

floor level. Since the two bands are very close to each other, 

Ids of both are the same. The ***** shows a bias of .boat 2.2 dB for an occupancy 
i^of 54 % relative to the low cha^l occupancy band. The variance also tncreases when 
t lupancy increase, In this p^culat case, *e variance is tripled. The »-*»-« * 

signals is available to the receiver. 

. v rttrt frhe noise floor level estimation 
Another embodiment of the noise noonc . 

„ - r- a wh i ch i s referred to herein as the penalized log-hkelibood 
method presented in Figure 6, which is Teierreuiu r 

r * method (601 602, 603, 604, 605, 606, 607) are the same 
method The first seven steps of the method (Wii, oiu, ou*, 

Tthe « seven steps .02, 1»3, I* KB. «. .«) of- DLCRtneurod of .he 
M bodiro*n, preyed ,a Pi^e 1. Thus. received sianal 60! to be analysed ~ed 
, to have been filtered to prevent ahastng of to freo.uet.cy bands after the sarupltng P»cess. A 
pluraluy of arraloeue-to-disita. converters «02 «**. the signal and stores a senes of 
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.aruules in . bufferreprescted by the vector r 603. The vector r is then windowed or 

L7e signs! P.* (noise and nowise signals, up to the fterbsn* « 

, anla input. * is — neee^v * cure that - - 

avoid any biased esbWion of the noise floor .eve!. ****** domam, the filte bar* 
oftenimplemented - aFastPounerTrausform ^because ofits — ^ 

integer log magnitude domain 606 ^ith a rounding operation or, the product of 20 tones the 

jt> va i ues in the vector x 607 are men uacu ^ ^ 

Itlpondin.to.hc^dB values. ^hU^^edisfhenused^a.^ 
!lTd serene* «. of abrade dB values m d« order w 609. The vector W rs Uren 
13 of the same long* and ~— the same _ as the vector s but with elemen* ^ 

equivalent to a decreasing sort operation on the rounded dB values. The preseM approach is 
much .ess complex in terms of number of operations recurred* perform the sorhn jta , a 

20 linear domain 610 with the transformjition 

v.-lO^ 0 

summed by groups of K610 srartmgvdth me largest values. This grouping reduces the 
dimensionof the problem from U. Another impact of the groups operation ,s *. 

tendency of the group sampie power, o be de.orrela.ed even wimthe ^""T; 
because me decreasing ordered samples are no, necessary from adjacent FET bras. The 
fonowmgolc^oUistnecompu^Uoooftop^^ 
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the expression 

PLLM(fc) = -L{k) + Cp{k) ,k = 0,1 , . . . , M - 1 , 

where L(k) is a quasi log-likelihood function of the M sorted group events, C is a constant 
a p«*y taction. In this ^ *e .uasi-log-liVcelihood Action „ ?1 ven 

by 



i(jt)=/Cln 



henotedtt^avanetyonog-.Utelihooiandqu^^^^^^ 

results. The function t « would be a true log..ikeli„e<,d iunenon 
variance of independent 1 «. —> O— variables,but the elements of a 

received srgnal are correlated »>d thus are no, indepentot nor arc tley necessartry zero » 
Onri. random Tb. penalty function,**) is added to the negative of the anas, 

w JMhood faction to form the function PLUW °v=r which the minim™ «11 be 
seLhed 614. The penalty fa«ta»PC*> - <*- — «° - *«-"*- 
« described in "Detection of Signals by Infcrrnaion Theoretic Criteria" M. Wax and T. 
Kailath, IEEE Transactions Acoustic, Speech, Sijoal Processing volume 33, numb,, 2, pages 
387-392, April 1985. The pena!ty function inthis embodiment, however, is different from 
those derived earlier intha, tie coefficients 4 are no. truly statistically independent, andfhat 
me intent is not to determine the fcD model order as in the reference cited above. For ttas 
embodiment, the penalty function is asecond order polynomial given by 
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p(k)=C 



OtherpCynomials or tart. a« p»ssible,b Ut *is secW orderpolynonu. - — foC 

signal environment a, the input. For to apphcatton, C - -2.6 » the preterr 

. , i that flic FFT coefficients art! distributed 

assumption to determine the above polynomial is that tram 

assump , . •„,„„,.,;„, If the additive noise time data samples 

according to a Gaussian probabinty density nmenon. It tneao .._,„. 
TJZc . Gaussian distribution the Unear combination of Oie FFT w,U tend to tt ausfonn 
£L> variable^ a probabiU* — on close to Oanssian v^ables ^a! 
m d large do of MK. The method is uierefore robust to non-idea, condtfons «*ere Uie 

noise might no, be trnlv Oaussisn. Th. scale ^r C sives a de.ee o Worn m 
fce optuni^onprocess. The output of Amotion PIXM(*> •. <*< f °'*< 

index of its minimum g KP defined as 

q m . = Wg m i n {- L ( k ) + Cp ^ ' 

The above expression means that q „ is the index at the minimum of the function 
PLLM(^ Tnenoisefloorestinaateper^FFTbiosistheav^ofthe M — q t p--l smallest 
numbers of the sorted vector values /,. The noise floor estimate 615 per ITT bin is then given 



by 



m. 1 

cr" = 



The noise floor level estimates can be averaged for the same frequency band to 
reduce the variance if the noise floor level is known to vary very slowly. 
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To „M the performance rf*. WW* -»-«»• *"*«'» ™^ 

ota-. Theperfom.ance in tennsoffcebiasand lb. vsrianceof thencse floor level 

of- .ethoo is presented . a funcuon of the .W oe.up.ncy and * a 

of the informal Si5 nal for . scenario of a speofic signal. A sampled FX 

si ™i in the 1 39 to 144 MHz band is used as the transmitted signs! in the smiulaaon 

Slgnalintne win dowed FFT flltet bank used a 

software The channel bandwidth is set to 30 kHz ana me * in 

to window. The transmitted sisnals all have the same average power a. the recover 

input 

H^mustratesmeno^^ 

signal. As the figure shows, the two curves are similar over the range of the index Figure 8 
presents the PLLM(*) function for a noise only, and a non-noise plus noise si^ai. It can be 
seeD that the index of the minimum of the PLLM(fc) function is indicative of the signal 
subspace pension. Since the remaning subspace includes noise only components, th«t the 
average is an estimate of the noise floor level. In the following, to reduce the amount of 
simulation results to present, the number of channels Af has been fixed to 6"4 and the 
of bins per channel K to 8. These numbers are typical of what is expected in practice. Other 
choices of number of channels and number of bins per channel are possible with Afferent 
performance results- Figure 9 illustrates the performances for the values of signaUo-no.se 
ratios of 10, 20, and 30 dB in 30 kHz with the true value of tbenoise floor level bong at 
0 dB Each column bar is the average of 1 000 independent estimates with an error b ar 
representing the standard deviation. Figure 9 shows that the estimate bias increases with the 
channel occupancy. The estimate is within 2 dB of the true value up to 84.375 % of channel 
occupancv 901 . Figure 1 0 presents real-time implementation results where 1000 estates of 
the noise floor level were computed using real-life signals. The signal of interest was over a 
4 5 MHz bandwidth sampled at 5.12 MHz with 5 kHz bin resolution. Two curves are 
presented: one where the centre frequency was at 848 MHz with a low occupancy of 
approximately 3 *. and one where the centre frequency is 880.05 MHz where the occupancy 
is about 54 % with an average signal power of 30 dB above the noise floor level. Smce the 
two bands are very close to each other, it is assumedthat the noise floor levels of both are the 
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same 



. xhe ^ - in— whea *e oeo»P=.cv 

go od considering <h* . . ^ Vi-W* of .he d** - 

^ „^ above arep^ded. * *. sate 

Ihe exempt y ™ rmir se numerous variations and 

, - . r - t t u irrotie of the invention. Ot course, nuni 

adaptations may be made to the aoove-aes 

depUs *- spirit of .he whi* is defined in Ale cfcnns. 

„ • ■ .h.^ecifically described embodiments can be earned 0« 

without departing fro™ the scope 01 ui«= m 
scope of the appended claims. 
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Window 


unset \oi>) 


Rectangular 


2 


Hanging 


2.75 


Bkcfctnan 


3 


Flat-top 


3.2 



Table 1 



